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Abstract 

We study twist expansions for parton fragmentation functions based on the definition of the 
twist as an invariant matrix element of a light-cone, bilocal operator. The results are then applied 
to a method which might be used to extract higher twist effects in the fragmentation sector using 
both e + e~ , and e~p collisions. We discuss how to apply the later measurements to experiments at 
the Jefferson National Acceleration Facility. 
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I. INTRODUCTION 



Duality in the nucleon structure function, W2{y,Q 2 ), was discovered before Quantum 
Chromodynamics (QCD) Thus the structure function, F 2 (Q 2 ^') = ^W^/m^y-, in the 
resonance region (W < 2GeV), is approximately the same as (duals) the deep inelastic 
region (W > 2GeV), when the functions are expressed in terms of the scaling variable, 
to' = 1 + W 2 /Q 2 . Here W is the final-state hadron mass. Moreover, the occurrence of 
duality appears to be local in the sense that it exists in each interval of u/ over the prominent 
nucleon resonances. 

An explanation of Bloom-Gilman duality was offered by de Rujula, Georgi, and Politzer 
in 1977 2j. Using the operator product expansion, they represented the scaling function as 
a sum of various twist operators, and then studied the contributions from individual twists 
(moments of the scaling functions) in the scaling variable, £ = 2x/(l + (1 + 4x 2 m 2 N / Q 2 ) 1 ^ 2 ) , 
where x = Q 2 jlm^v . They argue that the n th moment, M n (Q 2 ), of F 2 has the following 
twist expansion, 



M n(Q 2 ) = T,i-pyr) k ' lB nAQ 2 )- (i) 

k=l y 

Here, Mq is a mass scale (~ 400 — 500MeV) 2 , and B n k(Q 2 ) depends logarithmically on 
Q 2 , being roughly on the order of B n . According to eq. 1, there exists a region of n 
and Q 2 (n < jjp) where the higher twist contributions are neither large nor negligible, and 
where the dominant contribution to the moments comes from low-lying resonances. In 
this region for example, the moments defined by eq. 1 would not correspond to positive 
definite functions. Thus one might expect the structure functions to oscillate when large k 
is important, and this would lead to the appearance of local duality. A more recent study 
of parton-hadron duality by Ji and Unrau p| gives more quantitative estimations. 

There is no doubt that duality is a very interesting phenomenon, and it could allow one, 
under certain circumstances, to bridge the gap between perturbative predictions and experi- 
mentally observed quantities in non-perturbative regions [4]- 

However, duality, as expressed 

in the above analyses, only reflects properties of the parton structure functions. Appealing 
to the concept of factorization in the strong interaction, it is natural to ask whether there 
is duality in the parton fragmentation sector. Although in twist-two, perturbative calcula- 



tions, parton structure functions and parton fragmentation functions have some similarities, 
this does not indicate that they have similar higher twist expansions. 

Recently, an experiment at the Jefferson National Acceleration Facility, Jlab E00-108, 
studied duality in semi-inclusive ep reactions. Semi-inclusive reactions can reflect higher 
twist contributions in the parton fragmentation sector, so it is worthwhile to study these 
contributions in more detail in order to understand how one might extract these contributions 
in such reactions. 

We investigate these issues in this note. In Section 2, we study the twist expansions in the 
fragmentation sector, and investigate twist two perturbative calculations of semi-inclusive ep 
reactions. We also discuss a method to extract higher twist effects in these measurements. 
The paper is summarized in Section 3. 

II. TWIST EXPANSION IN THE PARTON FRAGMENTATION SECTOR 

In this section, we study the twist expansions of the parton fragmentation functions in 
e + e~ and ep semi-inclusive reactions. 

A. e + e~^h + X 

The kinematics for e + e~ — > h + X is illustrated in Fig. Two kinematic invariants, 
Q 2 and v — P ■ q, define the process. We consider scalar particles, h. The scattering 
matrix of the semi-inclusive reaction, e + e~ — ■> h + X, when calculated to lowest order in the 
electromagnetic interaction, is given by; 

e + (p e ') h(p h ) 



X(p x ) 
e(p e ) 

FIG. 1: Kinematics for e + e~ — > h + X. 




(hX\S QCD+ e m \e + e-) = {hX\Texp[ij d\{C^ + Cg CD )}\e + e-) 

Texp [i J d^Cg CD ] \e + e-) + 0(e 2 ) , (2) 

Here we have defined; 



W = T.Qi^M'h'Mp,)--, (3) 

where Q q is the charge of quark, q, in units of the proton charge, e. 

Using Feynman rules for QED, one then finds for the kinematics of Fig. Q 

M = ^v( Pe r, a e ,)^u(p e ,a e )(hX\J»(0)\0) . (4) 

Hence, the cross section can be written as; 

dcr ~ LuW^ - d ^ h ; (5) 



with; 



\{^v - Q 2 9^ ~ (Pe - Pe>)n(Pe ~ Pe')u} , and 

^ = ^E( 2 ^ 4 (^ + ^-'?)(o|^(o)|/ i x)(/ i x|^(o)|o) 



= hi d4 & iH J2(0\nO\hX)(hX\r(0)\0)^ (6) 

We have set the electron mass to zero. The sum over unobserved hadrons, X, cannot be com- 
plete because the state \hX) depends non-trivially on the observed hadron. Therefore one 
does not have £ x \hX)(hX\ = 1 and so Ex{0\MO\hX}(hX\J v (0)\0) ^ (0| J^)J u (0)\0). 
Thus e + e~ — > h + X is not controlled by the product of two operators, the operator product 
expansion does not apply, and no short distance analysis can be formulated. 

In a reference frame where the produced hadron, h, is fixed {(J, the hadron and photon 
momenta can be expanded as; 



1 , r~rr N ,. 1 



= —(v - V^ 2 " m»W + o (" + V" 2 " ™& 2 K , (7) 



where; 



Q 2 = (Pe + P 

Ph-q = v, 

%Ph ■ q 



2 

e'J , 



< Z 



q 



Writing; 



one finds in the Bjorken limit; 



P» = ^(1,0,0, 1), 

n" = —(1,0,0,-1). (8) 



f = np" + \n» + ^ , (9) 



lim q ■ £ = rtv — — . (10) 

Q 2 -*oo Z 



This result implies that as v — > oo, then r\ — > and A ~ z. Because £ 2 = £ 2 — £f — £ ±2 < 
£o ~~ £l> one finds in the Bjorken limit that £ M± — > 0. Therefore in this limit, light-like 
separation occurs, ~ 0, which dominates the integration region of eq. (6). 

Fragmentation is generally a non-perturbative process. We first consider the simplest 
quark fragmentation function represented diagrammatically in Fig. |21 More complicated 
fragmentation processes, such as coherent fragmentation of several quarks and gluons, do 
not contribute until order 1/Q 2 . 

In the case of Fig. EJ the diagram of e + e~ — > h + X shown in Fig. ^ can be re-drawn 
as in Fig. El 

The fermion field can be written as; 

1>V) = £ / j0^w [a{p > s)u{p > s)e " iP ' x + 6t(p ' s)v{p > s)eiP ' x] ' 

where a(p, s) {a^(p, s)) and b(p, s) (tf(p, s)) are annihilation (creation) operators for 
particles and anti-particles. Assuming that there is no final state interaction between Z 



FIG. 2: Quark fragmentation. 



h{Ph) 

= Y( P y) , 
- Z( Pz ) ' 

FIG. 3: Kinematics for e + e~ — * h + X for quark 6(6) with the Fig. |2] fragmentation, 
and hY, one can rewrite W^ v in eq. (6) as; Fig. |2] as: 

w»» = _L /" ^ e ;<K ^ (o| j^(OIM^ + ^)><M^ + ^)i^(o)|o) 

47r J yz{y+z=x} 

£ : {0|i.„K)|fty) 7 ^{0| T i fl (f)|Z>{Z|%(0)|0)7 i " A {fty|«)A(0)|0> : + 
yz{y+z=x} 

E = (0|^(OI^)7^(0|^(OI^ , )(^ / |^(0)|0) 7 ^(^Va(0)|0) :} 

Y',Z'{Y'+Z'=X} 

= rf 4 e[E(°i^(Oi^)7^(o|{^(e),^(o)}|o) 7 ^(/iy|^(o)|o) + 

E7^(0|^(OI^')(0|{^(O,^(0)}|0)(/ i F'|^(0)|0) 7 ^]. (12) 
y 

To obtain this result we have used the fact that; 




^2\Z)(Z\ = 1 and 

z 



^\Z')(Z'\ = 1 (13) 

Z' 



and due to the un-physical energy in the \Z) state, then; 
1 

= 5 4 (p z + g-^-p y )(0|^(0)|/ i F)7^(0|^(0)|Z)(Z|^(0)|0) 7 ^(^|^ A (0)|0) 



dV* £ (0|^(OI^)7^(0|^(0)|Z)(Z|^(O|0)7j A (^|^(0)|0) 
z,y{y+z=x} 



and : 



0; 



' ' d A ^< £ 7^(0|^(OI^)(0|^a(0)|Z)(Z|^(O|0)(^|^(0)|0)7^ 

Z,Y{Y+Z=X} 



= 6\p z + q-p h -p Y )Y a p(0\MmY)(0\M0)\Z)(Z\M0)\0)(hY\^ 

= 0. (14) 

In this expression, summation over color and flavor is assumed. 

Due to the strong interaction Lagrangian, Cqcdi ec L- (12) can be modified because of 
the higher order contributions in perturbative QCD calculations. These interactions cause 
vertex corrections, gluon polarization, etc, and their contributions to the cross section give 
logarithmic corrections. These contributions do not change the twists of the terms in the 
matrix elements. However, there are terms which represent quarks propagating in a gluon 
background in order to preserve color gauge invariance of the bilocal operator, W^ h ' . These 
terms are essential to generate higher twist corrections, and can be included by changing 
the singular function of the free field theory 

161 , ty(0, #))} = h MHom 2 ), to; 
MO, ^(0)} - ±- MZom 2 )V(exp[i / V^OD • (15) 
Substituting eq. (15) into eq. (12) and using the relationships; 

S^tpua = T-^ T 'y) 'pH pi 'vl 'a) 9p,p9ua 9p.u9vp 9ptu9up , 

(0\V(exp[i /V^VO)]|0> = J2(0\V(exp[i / ? d(»A»(())] \Z) ■ 
Jo g °° 

^•00 

(Z\V(exp[i rfC M ^(O)]|0), (16) 

J 

one has, 

1 r — r°° 1 

= -^J d^Yi^MOHexpl-i J dCA.iC))} \hX)—d p e{^)5(e) ■ 



(S^ J(T - te^ lal5 ) aX (hX\V(exp [i / d(»A,(()] )^ A (0)|0> + 



o 

oo 1 

1 r „ '^n/.2\ 



(0|^(0^(ea!P[i ^ rfC^(C)])l^)^[-^e(^ )5r)]- 



(S^T* - ie^w 5 ) Sf} (hX\P(exp [-i / dC4»(0] )^(0)|0> ] 

•/ 

n/b+n/b, (17) 



where Vu/biJ^h/l) ls t ne fragmentation function of a quark 6(6) fragmenting into hadron ft,. 

As shown in eg. (10), light-like separation ~ dominates the semi-inclusive process, 
e + e~ — ► h + X p. Therefore one can write; 

= An" + |", (18) 

where £ M includes contributions from rjp 11 and Both 77 and £ M± go to zero (l/A/Q 2 ") as 

g 2 -> oo 

With eq. (8) and n 2 = p 2 = 0, one can expand eq. (16) in terms of This gives, 

W^= J_ / d\ e - iX/z { 

x 

^[(oi^(An)p(expH / rf™-A(™)])|ftx)— a p 6(e )5(e 2 )- 

(5^ Tff - ie^ lal5 ) aX (hX\V(exp [i / drn ■ A{rn)\ )Va(0)|0) + 
(0|^(An)P(exp [z / rfrn • A(m)] )|ftA0— 3 p e(-eo)<5(£ 2 ) • 

/•oo _ 

(£*W" 7ff - ze^ 7(T75 ) 5/3 (/ i X|P(exp [-z / drn • A(rn)} )^(0)|0) ] 

JO 

+0(1)}- (19) 



The term in eq. (19) is the same as the definition of parton fragmentation given in 
ref.j]]]. Working in the light-cone gauge, n ■ A = 0, explicit reference to gluons disappears, 
and usingthe definition of the twist of an invariant matrix element of a light-cone bilocal 
operator p, ref. 8| studied the twist expansion of the possible terms for the production of 
scalar hadron or hadrons whose spins are not observable, 



and 



z J — e-^ z (0\Y^(0)\hX)(hX\^(Xn)\0) = ^{z)^ + / 4 (# 2 n"] , 

z / e -^/^o|^(0)|/iX)(/iX|^(An)|0) = AMeAz). (20) 



Here fi(z), e\, f± have twists 2, 3, 4, respectively; and M is a generic QCD mass scale. 

Since the twists for light-cone bilocal operators only represent the leading Q 2 dependence, 
it is possible for fi(z), t\ and/ 4 to include multiplicative factors of M 2 /Q 2 . Therefore eq. 
(20) can be rewritten as; 

z / ^e-^(0| 7 ^(0)|/.X)(/.X|^(An)|0)=4^[^(^)^n(^)+^(^)" + V4nW]; 
J 2tt n=Q VQ 

and ; 

z / -e~^(0\m\hX)(hXmXn)\0) = 4MY,(^) n e ln (z). (21) 
J ^ n= i vv 

One can see from eq. (17) that each £ factor always has a gauge-covariant derivative as a 
companion in the O (£-*-) terms; i.e. £ M -Dn- Using the light-cone, £ — > 1/y/Q 2 as Q — > oo and 
the twist analysis of bilocal operators p, one concludes that the C(£) terms have higher 
twists than their corresponding (£)° terms. For example, a general C(£) term of order n has 
the following twist expansion; 

f d\ 

Z I ' ' ' U 1 ^ ^ • • • D^M^)\hY)(hY\M0)\0) 

1 N 

VV 3=0 

(zAi* + fiV)M 2 ffi(z) + n ti n v M A f^~{{z) } + trace terms . (22) 

Based on the twist analysis from bilocal operators, the leading twists for fx+ii /jv+3 5 an d 
Jn+I are N + 2 + 2(JV - j), iV + 3 + 2(N - j), and N + 4 + 2(iV - j) respectively. 

Since all the trace terms have higher twists, their corresponding diagonal terms, eq. (22), 
can be written as; 

r rlX 

* J ^ e ~ iX/ ^i ■ ■ ■ C r ^(°l^ ■ ■ • D^M*n)\hY)(hY\M0)\0) 

N 00 yr 

~ J2p^ ■ .-p^n^ ■ ..n^[p^J2(^=) N+2 ^ N -^f^(z) + 

3=0 i=0 vv 

t^) N+2 ' N -' j)+l f N N+r^) + E(^f ) N+2{N - j)+i fS4>\z)] ■ (23) 

If there is final state interaction, there can be additional terms in W^ u . Ref. jsj has 
demonstrated the existence of such a terms, 

+ /• ^^ e iAA e <A.(iM-i/*)[ <0 |^(An)|7r(P)X)(7r(P)X|^(0)i5l(/xn)|0) . (24) 



In more complicated fragmentation processes, more quarks and gluons are usually involved 
in the operators, and these lead to larger dimensions. In any case, twist expansions can still 
be carried out based on a twist analysis scheme involving bilocal operators. 

In conclusion, one can carry out twist expansions for all bilocal terms in eq. (19), to 
obtain the twist expansion for W^ v \ 

W = E(^rTfn. (25) 

Twist two contributions to \M\ 2 are the major contributions as Q 2 —>■ oo or (^2) 3> n. 
As (j^) > n, higher twist contributions become important. Depending on the sign of 
the coefficients of (%^) n , the duality phenomenon in the fragmentation sector could occur. 
Although one cannot carry out non-perturbative calculations, the determination of the sign 
of f n could be measured in duality experiments. Such experiments are straightforward 
in e + e~ collisions, where one can compare the cross sections of hadron, h, production at 
fixed z for various Q 2 . For a given accuracy, e + e~ collisions at lower Q 2 require higher- 
twist contributions than those at larger Q 2 . However, it is meaningless to attempt a twist 
expansion in the non-perturbative region, n > (Q 2 /M 2 ). 

Thus the fragmentation functions can also be expanded in terms of contributions of 
operators of various twists, and the phenomenon of duality can appear if the signs of the 
coefficients of the higher twist contributions are not positive definite. Because complete 
non-perturbative calculations are not presently possible, a comparison of measurements of 
fragmentation at various Q 2 , having the same z, can provide information on the coefficients 
of higher twist contributions. 

Although we have studied only the production of scalar hadrons, a richer structure occurs 
if one includes polarization effects. 0,0]. 

B. The process e~p — > h + X 

The scattering matrix of the process e~ + p — > h + X is given by, 

W = (hX\Texp[z J d'i (£« + \ep) • (26) 

The lowest order of the electromagnetic contribution is shown in Fig. 0] 





VP 

FIG. 4: Semi-inclusive ep scattering. 

The scattering matrix can be written as 

W = $ e tfe)'f1>e(Pe)\(hX\J lt exp(i [ d^C£ CD )\p) + 0(e 2 ) , (27) 

q J 

where is the proton electromagnetic current, and depends on the proton structure. 
From eq. (26), one has, 

\W\ 2 = Mp' e )rMPe)Mp'e)rMPe) \ £<p|J>X) • (hX\J u \p) . (28) 

This is similar to e + e _ — > h + X, where Y^x \hX){hX\ ^ 1 due to the non-trivial de- 
pendence of \hX) on the observed hadron, h. One can work out a similar expansion as was 
done in the last subsection, to obtain the twist expansion of eq. (28). The major difference 
between eq. (28) and eq. (5) is that eq. (28) contains information on the parton structure 
functions in addition to information on parton fragmentation functions. 

The phenomenon of duality in parton structure functions has been established through 
inclusive ep scattering. In order to extract higher twist information in the parton fragmenta- 
tion sector using semi-inclusive ep scattering, one needs to separate the contributions from 
the parton fragmentation and structure sectors. If one allows factorization in the strong 
interaction, one should be able to separate the scale, Q 2 S , in the structure sector from the 
scale, Q 2 F , in the fragmentation sector. 

In general, it is more advantageous to obtain higher twist information in the fragmentation 
sector using e + e~ scattering. In order to extract information from ep scattering, one needs 



to subtract the contributions from the structure sector. This can be achieved by suitably 
choosing Q 2 S and Q 2 F in an experiment. One way to achieve this is to keep the Q 2 S sufficiently 
high so that the contributions from the parton structure sectors are mainly from twist 
two operators which can be confidently calculated. Then one can suitably choose the Q 2 F 
functions in the fragmentation sector to obtain information of higher twist contributions. 
Therefore to quantitatively study higher-twist contributions in the parton fragmentation 
sector, it is critical to identify Q 2 S and Q 2 F and separate the parton fragmentation sector 
from the structure sector. 

An experimental effort to detect higher twist effects in the fragmentation sector through 
semi-inclusive e~p — > h + X scattering has been recently undertaken at Jlab. We comment 
on this process in more detail here with the intent to clarify the experimental results. 

The object of such an experiment would be to determine higher twist effects by comparing 
measurements which include effects from all twists to the calculated twist two results. As 
discussed previously, in order to study higher twist effects in fragmentation sectors through 
semi-inclusive e~~p —>■ h + X scattering, it is important to subtract the contributions from the 
parton structure functions. This is feasible by choosing high Q 2 S in the structure sector so 
that these contributions are predominantly from twist two operators which can be calculated 
perturbatively by using the well-known information on twist two parton structure functions. 
One must have measurements in a sufficiently high region of Q 2 S , with various Q 2 F , in order to 
compare the measurements with calculations of twist two contributions in the same regions. 
The differences between the measurements and the calculations provide information of higher 
twist effects in the fragmentation sector. Therefore we need to locate the regions with high 
Q 2 S and various Q% to carry out the twist-two calculations. 

If Q 2 S is large enough so that the contributions from the parton structure sectors are twist 
two, the scattering process e~p — > h + X is illustrated in Fig. El 

Before presenting a calculation of the process shown in Fig. we would like to make several 
remarks. 

1. In a momentum infinitive reference frame where \p\ 3> m p , parton momentum and 
energy are usually expressed as a fraction x of the momentum (energy) of the parent 
proton momentum (energy), where x = m the CMS or Lab frame, a Jlab beam 
energy (E e ~ 5.5 GeV) is too low to the momentum fraction since 

the beam momentum p e in both CMS and Lab does not satisfy \p e \ ^> m p . In our 



FIG. 5: Semi-inclusive ep scattering Kinematics for e + e — * h + X. 

approach in this section, we work in the target rest frame. While the target mass, m p , 
is non-negligible when compared to the beam energy, the quark masses are. Therefore, 
the quark energy and momentum in the target rest frame can be expressed as; 

\p g \ = xm p . (29) 

In a finite momentum reference frame, parton transverse momentum cannot be ig- 
nored and the parton structure function f q should depend on x, Q 2 , and p^, i.e. 
fq = fq(x, Q 2 ,Pq), where x = E q / Etarg- This structure function should be the same as 
the one measured in an infinite momentum reference frame, if one boosts the reference 
frame to a very high momentum so that pT is negligible. In a target-rest reference 
frame, the probability for a parton to have a fraction x of the parent proton's en- 
ergy should be the same as the corresponding parton structure function measured in 
momentum- infinitive reference frame because x = E q /E targ . However, the direction of 
the momentum is completely random and satisfies Y^,qPq = 0. 

2. In e + e~ scattering, the Q 2 F is defined as Q 2 F = (p e + + p e -) 2 and is equal to 4E 2 eam in 
the CMS. If one concentrates on the fragmentation functions given by Fig. the 
energy of the produced quarks, E q , is equal to E beam in the CMS. Therefore Q 2 F can 
be written as; 

Q 2 F = AE 2 q . (30) 

The fraction of energy carried by the produced hadron from the parent parton is 
z = 2p h ■ q/q 2 , which is equal to E h /E q in the CMS. 



In semi-inclusive e p scattering, one usually defines (see Fig. 

Qp = ~Qs ' 
Pp-Ph 



101 



Z = 



pp ■ q 



(31) 

lab. frame 

In order to adopt the fragmentation functions obtained from e + e~ scattering in the 
e~p scattering process, i.e. 

= VtjfiQ'lz'), 

one should use the same definitions used in e + e~ scattering. Therefore, we correspond- 
ingly redefine the Q 2 F and z in e~p scattering as 

Pp -Ph 
z = 

Pp ■ Pq> 

-(¥) ■ 

\ 1 / lab. frame 

2 ^PP-Pq'? 
VF — 2 
PP 

= UeI) . (32) 

V 1 J lab. frame v ' 

Compared to Qp and z in e~e + scattering, we believe that these definitions are rea- 
sonable and the new Qp is at least proportional to Q\ in e~e + . 

3. The Q 2 S in ep scattering is defined as usual as; 

Qs = -fs = -(Pe-Pei? ^ 2E e E e ,(l - cos0 ee ,) . (33) 
The differential cross section of the process e~p — > h + X is given by, 

Me-p^h+x) = £/^w§ |ge _ g ; 2ge2gg ^i^- 

where p e , p e r, p q , and p g / are four momenta of the incoming electron, scattered electron, 
quark within a target proton, and scattered quark, respectively. Also f q (x) is the structure 



function of quark q in a target proton with x fraction of the proton energy, D q (Q F , z) is the 



fragmentation function of quark q fragmenting to hadron h with z fraction of the quark q 
energy, Q q is the solid angle of quark q within the target in the lab frame, and the 5-function 
reflects energy-momentum conservations among p e , p e i, p q , and p q i. 

We define the electron beam direction as the positive z direction and let; 



O 2 

W 1 = E e E q , (1 - cos 9 q ,) + ^(E e -E q , cos . 

zhp 



'q ) > 



W 2 = WrE^l - cos 9 q/ ) + ^-E 2 sin 2 9 q , cos 2 ( 



W 3 = W 2 + 

Eq. (34) can be rewritten as; 

da{eTp — > h + X) 
dQ 2 s dQ 2 F dz 



where 



Q 



s ^ 2 El sin 2 9 a , cos 2 ( 



2E e ' 9 



(16n) 3 m 2 P ^ 



x 2 (l — cos 6^) 



e /d cos 9 ' q 'd(f) q t 



\M\ 2 D h (Q 2 F ,z) 



v e . - v, 



q\E 2 E q 



■fq(x) ■ 



E P i 



W 2 + JW2 - W 3 E 2 ,(1 - cqs6 ,] 



E z g ,(l-cos9 q ,] 



Wo 



W 2 - W 3 E 2 (1- cos 9 q i) 
E 2 (l-cos9 q ,) 2 



if cos( 



if cos( 



-<M > 



- M < o 



x 



cos 9 n 



sin 9„ 



COS( 



sm i 



E e i + -Eg' — -E e 

irtp 

Epi COS # e / + E q i COS 6^/ — i? e 

xmp 



E 2 sin 2 # e / + E 2 , sin 2 6*5/ + 2E e >E q i sin 6* e / sin ? / cos(0 e / — 4> q >) 

x 2 m 2 P 

E e i sin 9 e i cos e / + .Ey sin 9 q i cos c/y 



im p sin #5 
£v sin 9 e i sin e / + £y sin 9 q i sin g / 



(35) 



(36) 



(37) 



xmp sin g 

Using t he p arton structure functions from ref. jl(| and the pion fragmentation functions 



;ne pt 



from ref. the pion production cross sections versus z at Q 2 S = 10.0 GeV 2 and £y = 

5.0 GeV and 4.0 GeV^ are plotted in Fig. |H1 In general, the production cross section is small 
and the choices of E q i that satisfy all conditions are limited. This can be seen from the very 



small cross section at E q > = 4.0 GeV in Fig. |H]when Q 2 S = 10.0 GeV 2 . This problem can 
be avoided by lowering the beam energy, which leads to lower Q 2 S and more choices of E q i. 
However, lower Q 2 S brings more contributions from higher twist operators in the structure 
sector. When the magnitude of Q 2 S is sufficient that the parton structure functions are safely 
represented by twist 2 operators, one can compare the hadron production data at various 
Eg/ to obtain information on higher twist contributions in the fragmentation sector. 




E; = 4.0 GeV 
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Z 

FIG. 6: Calculations of da / dQ 2 s dQ 2 F dz versus z for Q 2 S = 10.0 GeV 2 and E q > = 5.0 GeV and 
4.0 GeV. 



The expressions given in eq. (34) are given for the purpose of calculating the hadron 
production cross section versus z at Q 2 S and E q i. In experiments, the quantities that can be 
directly measured are E e >, 8 e >, 4> e >, Eh, and approximately 9 q > and <f) q >, if one takes 9 q > ~ 9h 
and (f> q > ~ (ph- The quantities that can not be directly measured are x, 9 q , (fi q , E q >, and z. 
One can determine x, 9 q , <f) q , and E q i based on energy and momentum conservation and the 
measurements of E e ', 9 e ', <fi e >, 9 q /, and <fi q > with the following equations, 

Q 2 S = 2E e E e/ (l-cos9 e ,), 

E t _ E e E e i(l - cos9 e ,) 

q E e i{l — cos9 e > q r) — E e {l — cos9 q i) ' 

E e > + E q i — E e 

x = , 

mp 

E e i COS 9 e i + E q i COS 9 q i — E e 

cos 9 e = , 



m P 



sin# e 



{E 2 e , sin 2 6 e i + E 2 , sin 2 6 q > + 2E e ,E g , sin 6 e > sin 6> 9 / cos(0 e / - 4> q ') 1/2 



cos e 



E e , sin # e / cos e / + E q , sin ? / cos c/y 



sin e 



mp sin 6 e cos0 e 
-E e / sin # e / sin <p e , + E q , sin 9 q , sin </y 



cos#, 



mp sin^e sin e 
cos 6* e / cos 6 q i + sin # e / sin 6 q i cos(0 e ' — <fi q >) . 



(38) 



Also z can be determined through the above determination of E q , and the measurements of 



With eqs. (38) and (39), one can determine Q 2 S , Q 2 F , and z of the hadron production 
processes and then obtain the cross section versus z. Comparing the cross sections at various 
Q 2 F but with the same Q 2 S and z, one can deduce the higher twist contributions from the 
fragmentation sector. 

III. SUMMARY 

We have studied possible higher twist QCD contributions in the fragmentation sector. In 
summary we find the following. 

1. The contributions to e~e + — > h + X cross section can be expanded as contributions 
of operators with various twists. If Q 2 F — > oo, twist-two operators are the dominant 
contribution. As one decrease Q 2 Fl higher-twist contributions become more and more 
important. If the sign of the coefficients in higher- twist contributions are not positive 
definite, oscillation similar to duality in e~p inclusive scattering could appear. As Q 2 F 
becomes low, QCD becomes completely non-perturbative, and the twist expansion 
becomes meaningless. 

2. A study of higher-twist contributions can be carried out through e + e~ — * h + X by 
comparing measurements of cross sections versus z at various Q 2 F . In order to carry 
out the same study using e~p — > h + X, one needs to make sure that the higher twist 
contributions from the parton structure sector are well determined, so that one can 
isolate the contributions from the fragmentation sector. One way to achieve this is 




(39) 



to choose large enough Q 2 S so that the contributions to the cross section from the 
structure sectors are dominantly from twist-two operators, which can be calculated 
with well measured structure functions and perturbation QCD. To achieve higher Q 2 S 
but different Q 2 F , one can suitably choose the detector angles carefully in e~p — > h + X 
experiments for various beam energies. 

3. We have not attempted to study the evolution of the ra,tio,a e - p ^ h x(Qs) / cr tot (Q 2 s ) as 
a function of z. However we suggest that this ratio contains information of the higher 
twist effects in the fragmentation sector and should be studied in the future. 

4. To adopt the fragmentation functions obtained in e~e + — * h + X in the process of 
e~p — > h + X, redefinitions of Q 2 F and z are necessary. We currently use the definitions 
in eq. (31). 

More studies on this subject, especially considering the case of polarization, are necessary. 
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